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We calculate the absorption spectra of a semiconductor microcavity into which a non-equilibrium 
exciton population has been pumped. We predict strong peaks in the spectrum corresponding to 
collective modes analogous to the Cooper modes in superconductors and fermionic atomic gases. 
These modes can become unstable, leading to the formation of off-equilibrium quantum condensates. 
We calculate a phase diagram for condensation, and show that the dominant instabilities can be 
at a finite momentum. Thus we predict the formation of inhomogeneous condensates, similar to 
Fulde-Ferrel-Larkin-Ovchinnikov states. 



I. INTRODUCTION 

The appearance of order at an equilibrium phase tran- 
sition is a central concept in many areas of physics, from 
condensed matter to the physics of the early universe. 
Recently there has been considerable interest in the more 
general problem of ordering far from thermal equiUbrium, 
motivated by the possibility of quantum quench exper- 
iments in cold atomic gases.^^ In a quench the parame- 
ters of the system are rapidly switched from a disordered 
to an ordered phase, and the disordered state forms the 
initial conditions for a dynamics with the new param- 
eters. An interesting regime is that of coherent relax- 
ationless dynamics, which can lead to the formation of 
non-equilibrium order including crystallization,'^ conden- 
sation, and ferromagnetism.^ 

Among condensed matter systems, semiconductor mi- 
crocavities are promising candidates for studying such 
quench dynamics. The non-equilibrium dynamics of 
microcavities has attracted considerable interest both 
experimental!}''^'^ and theoretically,"^'''^ with recent ex- 
periments demonstrating regimes where the low-energy 
quasiparticles, polaritons, form a condensate. More re- 
cently, an experiment has been proposed to implement a 
quantum quench,!^ by rapidly preparing a microcavity in 
a non-condensed initial state. The coherent dynamics of 
this non-condensed state is predicted to lead to a form of 
non-equilibrium condensation, similar to that predicted 
in a quenched Fermi gas.'^^H^^ We show here that, as in 
the Fermi gas, such condensation is due to the appear- 
ance of a new collective mode. Moreover, we show that 
in the microcavity the dominant instability occurs at a 
finite wavevector. Thus we predict that microcavities 
could be used to realize inhomogeneous condensates,'!^! 
i.e., those characterized by a spatially varying phase. 
These condensates are similar, in this essential respect, to 
those predicted by Fulde, Ferrel, Larkin and Ovchinnikov 
(FFLO) in unbalanced Fermi systems.'iSl 

In this paper, we first calculate the optical spectra of 
a microcavity a short time after it has been prepared 
in a non-condensed state, i.e., immediately after the 
"quench". We find that the collective mode responsible 
for condensation is directly observable in these spectra. 
We use this analysis to calculate a phase diagram for the 



non-equilibrium condensation, and show that the con- 
densation generally occurs at a finite momentum. While 
we focus on a microcavity containing quantum dots, our 
analysis is based on the Maxwell-Bloch equations. These 
describe a wide variety of coupled light-matter systems, 
implying a broad relevance of our work. 

The remainder of this paper is structured as follows. 
In Sec. [n] we briefly review the proposed quench exper- 
iment, and outline our model. In Sec. |III| we present 
absorption spectra of the system. In Sec. |IV| we discuss 
the phase diagram and the possibility of finite momentum 
condensation, and in Sec.|V]we discuss the connections to 
FFLO and the role of nonlinear terms. Section |VT] sum- 
marizes our conclusions. Finally, the appendix contains 
a brief treatment of the preparation of non-condensed 
initial states by optical pumping. 



II. MODEL 

We consider an experiment, proposed in Ref. I15[ on 
a set of localized exciton states in a planar semiconduc- 
tor microcavity. Such excitonic states could be realized 
in practice using either highly disordered quantum wells 
(where excitons are localized by disorder) or quantum 
dots. The proposed experiment involves two stages which 
are separated in time and can be regarded as indepen- 
dent. In the first stage, the localized states are driven 
by a chirped laser pulse. This pulse creates an energy- 
dependent population in the inhomogeneously broadened 
exciton line by adiabatic rapid passage (ARP). For cer- 
tain populations a second stage may then occur, where 
the population evolves into a non-equilibrium condensate 
due to the photon-mediated interactions between the ex- 
citons. 

As in Ref. [K] we describe the system using a general- 
ization of the Dicke model.'22l The localized exciton states 
are treated as two-level systems, with the standard dipole 
coupling to the electromagnetic field. The state localized 
at site i is specified by the Bloch vector ai — (o'i), where 
the inversion = 1(— 1) for an occupied (unoccupied) 
state, and a~ is the exciton annihilation operator. An- 
gle brackets () denote expectation values in the quantum 
state of the system. 
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We consider time scales short compared with the ex- 
citon lifetime, which is at least 100 and treat the 

electromagnetic field using a mean-field approximation. 
In this approximation the photon creation and annihila- 
tion operators are replaced with their expectation values, 
and hence become c-numbers. The resulting equations of 
motion are linear in the remaining operators, so that we 
may take their expectation values without further ap- 
proximation. The resulting dynamics obeys the general- 
ized Maxwell-Bloch equations: 



iipk = '^ki/'k + 9 



iPk = EPk - .9 X! ^k-k'V'k' , 

k' 

iDu = 2.g^ (P^^-kV-k' - Pk'+kV-kO 



(1) 
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Here V'k is the complex amplitude of the microcavity 
mode with in-plane wavevector k and energy Wk = !)• 
It is related to the expectation value of the photon anni- 
hilation operator by ■(/'k = (V'k) / v'^, where N is the total 
number of localized states. This normalization is conve- 
nient when dealing with condensation, since macroscopic 
occupation corresponds to a finite 'i/'k in the thermody- 
namic limit iV — > oo. We allow for the finite lifetime of 
the photon modes by taking 3cJk = —7- /k is introduced 
to allow us to calculate the linear response. F is an ex- 
ternally applied pump field, with wavevector p, that is 
used to create the non-equilibrium population. 

The coupling g = giy/n in Eqs. (1 3) is related to the 
dipole couphng strengths of the localized states, gi, and 
their area density n. To simplify the notation we have 
taken gi to be the same for all states; the extension to a 
distribution is straightforward. In the dipole gauge 
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where d is the matrix element of the dipole operator ef 
between the zero-exciton and one-exciton states, and Ei 
their energy difference, w is the effective width of the 
cavity, which arises from the normalization of the cavity 
mode functions. 

Pk{E) is the collective polarization of the ensemble at 
wavevector k due to states with energy in a small interval 
near E, 



Pk(i?)5i?=^^'(a-)e 



(5) 



where the prime indicates that the sum runs over states 
with exciton energies between E and SE. D]f^{E) is the 
collective inversion, defined in a similar way with re- 
placing a~. 

In the following we shall be concerned with large N, 
and wavevectors which are small compared with the in- 
verse spacing of the localized states. In these limits we 



may approximate sums over dot positions, such as those 
in Eq. ([5|, by 
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Here is the distribution of localized states in en- 

ergy, normalized to one. Thus Nv{E)5E is the number 
of terms in the primed sum, Eq. For k ^ the pha- 
sor sum is a 2D random walk, producing the 0{l/^/N) 
corrections.'^ 

The approximation of Eq. ([6]) corresponds to replac- 
ing the response of the disordered dielectric with its 
homogeneous average response, so that the wavevector 
is well-defined. This is similar to the linear disper- 
sion model which has been extensively used for inor- 
ganic microcavities.^"^ The 0{l/^fN) corrections describe 
Rayleigh scattering from density fiuctuations in the di- 
electric. They are generally small corrections because 
the scatterers are dense, so that on long wavelengths the 
medium appears homogeneous. The corrections ca n be- 
come important at very small or large wavevectors 12il2ll 
for modes whose group velocity becomes very small. In 
this case a long lifetime is required for the wavevector to 
be well-defined, so that even weak scattering or absorp- 
tion destroys the quasi-propagating modes. Here, how- 
ever, we are concerned with modes that have a significant 
dispersion due to their photon component. Furthermore, 
the lifetime of these modes is massively enhanced by reso- 
nant gain from the populated excitons. Thus the leading 
approximation of Eq. (|6| will capture the physics at the 
experimentally relevant wavevectors. 



III. ABSORPTION SPECTRA 

In the experiment is a chirped pulse, which creates a 
non-equilibrium population of excitons using ARP. In the 
appendix we demonstrate this explicitly, using a model 
pulse, Eq. (A5|, for which analytical solutions to the dy- 



namics exist. Following the pulse the exciton states are 
populated with a distribution given by Eq. ( A6 1 and the 
fields and polarizations are negligible ^k ~ 0, Pk ~ 0. 

To establish the optical properties of the microcav- 
ity immediately after the pump pulse we find the re- 
sponse to a weak probe /k. The susceptibility can 
then be found from the induced electromagnetic field 
(5V'k = X)k' / Xk,k'(^ — t')fk'{t') dt'. If the system is sta- 
ble then ■(/'k and Pk are small (of order /k) for all times, 
whereas if it is unstable they are only small soon after the 
pump pulse. In both regimes we may neglect terms above 
first order in i/'k and Pk- Eq. (|3| then gives I)k = 0, 
so the non-equilibrium population is constant. Fourier 
transforming the linearized Eqs. ([!]) and ([2| gives 

cjJVk = Wk^V'k " / ^^^^ dE + fi^. (7) 



The pumping populates the states independently of 
their position, within the pump spot. Thus the sum in 
-Dk-k' is strongly peaked near the forward scattering di- 
rection k — k' = 0, as discussed above [Eq. Neglect- 
ing the smaller off-diagonal scattering terms we obtain a 
diagonal response function 
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The absorption coe fficien t of the microcavity follows 
from the susceptibilitj^^EIl 



A{io) = —2 lim Imx(a; 

e— >-o 



ie), 



(9) 



where the infinitesimal e appears due to causality, and 
can be physically understood as a small damping con- 
stant for the excitons. The sign is such that A{uj) > 
corresponds to absorption of energy by the system. Thus 
from Eq. ([8| we obtain 



A{uj) = 2 



= 2 



7 - 5^7rDo(a;) 



(10) 



When the dots are unoccupied -Do(w) < and the empty 
exciton states contribute to absorption. For energies 
where there are occupied exciton states Do{uj) > 0, de- 
scribing gain due to the population. 

In general the response, Eq. (10), peaks near the ze- 



roes of G{u}), which are at the energies of the normal 
modes. These energies differ from the energy of the cav- 
ity resonance due to the coupling to the exciton states. 
For an unpopulated state the condition G{ll!) = recov- 
ers the usual polariton energies of the Lorentz oscillator 
model,'2Sl but in general the spectrum differs due to the 
presence of the non-equilibrium population. The modes 
have a lifetime determined by the second factor in the de- 
nominator, with contributions from the cavity losses and 
the resonant mixing with the band of exciton states. As 
expected it is damping which controls the overall strength 
of the absorption, so that the damping factor H{uj) also 
appears in the numerator. 

The only dependence of the spectra, Eqs. ^ and (10 1, 
on the wavevector is in the energy of the cavity mode, 
ujc = ^i^k) « ciJo + |kp/(2m). We therefore show results 
as functions of cUc, which corresponds experimentally to 
both the incident probe angle and the cavity width. 

Figure [T] illustrates absorption spectra obtained from 
Eq. ( |10| for both a pumped and unpumped exciton 
line. These spectra are valid at all times if conden- 
sation does not occur (see later), but only soon after 
the pump pulse if it does. We have taken a Gaussian 
model for the inhomogeneously-broadened exciton line, 
with standard deviation a, and measure energies relative 
to the center of the line. We choose the duration r of 
the pump pulse as our unit of energy, and have taken 




(a)Populated 




(b)Unpopulated 



FIG. 



1: Absorption spectra [equation (lOl] in arbitrary units, 
for the parameters described in the text. Top panel: spectra 
with the exciton population created by the pump pulse. Bot- 
tom panel: spectra of the unpopulated microcavity. Each 
curve corresponds to a different value of the cavity mode 
energy lJc, vertically offset as indicated by the right-hand 
axes. Note the peak developing in the populated system 
at tjprobc ~ — 13/r, indicating the presence of a collective 
mode. There are regions of probe gain {A{u}) < 0) between 
—31 < ujprobeT < "13, whcrc gain from the populated exci- 
tons overcomes the cavity losses. Dotted lines indicate spectra 
with unstable normal modes. Instabilities occur when a nor- 
mal mode lies in the region of gain and in these cases the ab- 
sorption spectra have negative peaks, which are hidden from 
view in this figure. 



g — 13/t, 7 ~ 1.5/t, cr — 15/t. These parameters, with 
r = 3 ps, are reasonable for a microcavity containing 
strongly disordered quantum wells.'^s! As discussed in the 
appendix the pump creates a population equivalent to a 
Fermi function with temperature l/^TTkgT), and Fermi 
energy fi dictated by the chirp and center frequency of 
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the pulse; we choose a pulse for which ^ = — 12.5/r. 

The lower panel of Fig. [T] shows the expected result 
for an unpopulated microcavity. There is a pronounced 
peak in the absorption at the cavity mode energy, which 
broadens as the cavity mode is tuned through the exci- 
tons. There is some suggestion of an anticrossing near 
resonance, i.e., a polariton splitting, but since the inho- 
mogeneous broadening is relatively large compared with 
the coupling this is a weak effect. The top panel shows 
that the population dramatically changes the absorption 
spectrum. For these parameters it leads to a range of 
probe frequencies for which 7 < 5^7r_Do(w), and the ab- 
sorption coefficient, Eq. 



(10) 



becomes negative. This 
occurs when the gain from the populated exciton states 
exceeds the losses, so that there is a net gain for the probe 
beam. Moreover, we see a pronounced additional peak 
in the absorption spectrum, which first appears near the 
upper edge of the population as the cavity mode energy 
is decreased. As the cavity energy is decreased still fur- 
ther this peak moves down through the region of gain, 
before the spectrum finally reverts to one dominated by 
the unperturbed cavity mode. 

This additional peak in the absorption spectrum is 
analogous to the Cooper pairing mode in a superconduc- 
tor or Fermi gas, that gives rise to the Cooper instability. 
The analogy can be seen by noting that the normal-mode 
condition G{uj) = contained in Eq. ( 10 ) is the Cooper 



equation, as discussed for this system in Ref. The 
non-equilibrium exciton population corresponds to the 
Fermi distribution, while the photon-mediated interac- 
tion between excitons corresponds to the pairing interac- 
tion between the electrons. As in a superconductor the 
sharp step in the population leads to collective modes 
generated by the pairing interaction. Fig. [l] shows that, 
for reasonable parameters, these collective modes give 
rise to strong features in the spectra. 

It is interesting to compare the spectra of Fig. [T] with 
the predictions for an equilibrium condensate in the same 
model.'23 In that case the condensation opens a gap in 
the single particle spectrum, which is the analog of the 
Cooper gap of the superconductor. Inside this gap is a 
collective mode, which is the analog of the Cooper mode 
or phase mode of the superconductor. The features visi- 
ble in Fig.[l]arise from the non-equilibrium generalization 
of the collective mode (which is a different spectral fea- 
ture than the gap). It is clear from Fig. [I] that it is the 
collective mode which dominates the spectrum. Thus, 
although the single particle features may be affected by 
condensation, this would have little effect in practice. It 
may be possible to isolate the single particle spectrum in 
a Rayleigh scattering experiment, as has been proposed 
for equilibrium condensates .'221 



IV. PHASE DIAGRAM 

The normal modes of the system, with frequencies de- 
termined by G{uj) = 0, have decay rates H{uj). If a 



100 
80 
60 
40 
3" 20 


-20 
- 40 



Normal 



Condensed 








10 



20 



30 



40 



50 



60 



FIG. 2: Non-equilibrium phase diagram for a populated mi- 
crocavity with a single photon mode as a function of cou- 
pling strength g and cavity detuning ujc with cavity damping 
7T = 1, 15, 30, 45. Arrow indicates curves of increasing 7. 
Dotted line indicates the location of the equilibrium phase 
boundary,^^ with temperature and chemical potential corre- 
sponding to that of the pumped population [T = l/(fcs7rr), 
11 = -12.5/r]. 



normal mode frequency lies in the H[lj) < region pro- 
duced by the non-equilibrium population it will be unsta- 
ble, growing exponentially to give a state with a highly 
populated mode, i.e., a condensate. The condition for 
the onset of such an instability gives a non-equilibrium 
phase diagram, which is shown for our chosen parameters 
in Figs. [2] and [a] 

Fig. [2] shows the phase diagram assuming that only a 
single cavity mode, of energy ujf,, is relevant. The dot- 
ted line shows the phase boundary for equilibrium con- 
densation in the same model, with a temperature and 
chemical potential corresponding to the pumped popula- 
tion. We see that one sheet of the non-equilibrium phase 
boundary extends the equilibrium result to allow for the 
cavity damping. Whereas in equilibrium the presence of 
the collective mode is sufficient to create an instability, 
in the open system condensation only occurs if the gain 
at the energy of the collective mode overcomes the cav- 
ity loss. Thus the collective mode can exist even in the 
normal state (see Fig. [T]), and the damping pushes the 
transition to larger couplings. In addition, we see that 
there is a lower limit on Wc in Fig. [2] This lower threshold 
is a purely dynamical effect, not present in the equilib- 
rium case. Below it there is a bosonic collective mode 
at an energy well below that of the populated states (see 
Fig. [1]). Although this mode would be occupied in equi- 
librium it is far out of resonance with the excitons. As 
a result, it is not occupied dynamically, and the uncon- 
densed state is metastable. A similar metastable region 
has been predicted in quenched atomic gases.'i^ 

Figs. [1] and [2] show that, for a given g and 7, the 
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FIG. 3: Non-equilibrium phase diagram for a populated mi- 
crocavity with a continuum of in-plane photon modes as a 
function of the coupling constant g and detuning at k = 0, 
(jjo- The cavity damping 7 = 1.5/r. In the inhomogeneous 
region modes with k 7^ have the highest growth exponent. 
The vertical boundary is a result of the gain-loss criterion, 
H{lj) > 



condensation instability occurs over a range of ujc- In 
a microcavity different values of lOc = '^o + |kp/2m 
correspond to either changing the cavity width, which 
varies the detuning loq, or considering modes at a differ- 
ent wavevector k. As such, a range of unstable oJc implies 
that for a fixed cavity detuning there can be instabilities 
at many wavevectors, with different growth exponents 
At short times after the population has been 
created the mode with the highest growth exponent will 
dominate. Figure |4] shows that as loq is lowered this dom- 
inant mode occurs at k 7^ 0, implying a condensate with 
finite momentum, and a spatially inhomogeneous order 
parameter. Thus the full phase diagram, allowing for the 
continuum of in-plane modes, takes the form shown in 
Fig.|3] 

The phase diagram of Fig. [3] can be understood phys- 
ically by noting that the condensation is a result of the 
exciton-photon interactions. If the cavity mode is de- 
tuned well below the excitons then the quasiparticles at 
k = are essentially photons, uncoupled from the exci- 
tons. Thus the condensation shifts to the higher momen- 
tum states, where the photons and excitons are closer to 
resonance, and there are strong coupling effects. 



V. DISCUSSION 

There has been extensive theoretical work on states 
with finite momentum Cooper pairing in the context of 
equilibrium superconductors, atomic gases, and quark 
matter.'iSl These FFLO states, which have been sought 
in a wide range of systems, may be the ground states 




FIG. 4; Growth exponents —H{(jj) of the unstable modes with 
wavevector modulus A; in a populated cavity with detuning ujq, 
coupling strength g = 13/r, and decay constant 7 = 1.5/r. 
The wavevector of the most unstable mode increases as ujq is 
reduced. For cjc ^ — SO/r the k = mode is stable and there 
are only finite momentum instabilities. 



where there is an imbalance in the populations of the 
two pairing species. However, they involve increasing 
the kinetic energy in order to gain pairing energy, and in 
practice this restricts them to small regions of parameter 
space. Here, however, the state achieved is determined 
by the Cooper equation and a gain-loss criterion, with 
the energetics playing a subsidiary role. Thus, as indi- 
cated by Fig. |3j condensation at a finite momentum may 
be achieved without fine tuning of parameters. 

The connection to FFLO may be made more explicit 
using a representation for the exciton operator a~ = 
Ci.-fCi,^ where c is a fermionic annihilation operator. In 
the simplest case of a plane wave condensate at wavevec- 
tor k the mean- field order parameter is P^, the macro- 
scopic component of the exciton polarization. This be- 
comes 

i 

p,q i 

~ (l/iV)^(cp,tc_p_k4). (11) 
p 

Thus we see that the condensate can be formally rep- 
resented as a coherent state of fermions, pairing with a 
finite total momentum. While in this respect the state 
is similar to FFLO, there are other differences. For ex- 
ample, in Eq. (11) the relative wavefunction of the pair 



is independent of momentum, and the pairing is entirely 
local. In a general FFLO state there is a momentum 
dependent pairing function, describing Cooper pairs of 
finite size. 

Because the growth exponent depends only on |k| the 
condensate emission at short times will cover a circle 
of in-plane wavevectors, giving a cone of emitted light. 
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However at later times the nonlinear terms neglected in 
Eq. ([7]) will break the degeneracy, selecting a spatial 
form for the condensate. In equilibrium such interac- 
tions favor condensate structures consisting of a pair of 
antipodal wavevectors (k, — k), or more complex struc- 
tures such as face-centered cubes.'^ Here the nonlinear- 
ity corresponds to the depletion of the exciton popula- 
tion by the growth of the condensate. This will reduce 
the gairli^l for collective modes of similar energies, sug- 
gesting that a single plane wave (Fulde-Ferrel) state may 
be favored. Although these nonlinearities determine a 
particular form for the condensate it is unlikely they will 
lead to a homogeneous state, so we do not treat them in 
detail here. 

It is interesting to note that finite-momentum polari- 
ton condensates have been observed,** ^ ^ though in a dif- 
ferent experimental protocol to that considered here. In 
these cases there is continuous pumping and relaxation, 
and a spatial structure imposed by a pump and trap. 
The mechanisms leading to this finite momentum con- 
densate have yet to be established, and are likely dif- 
ferent from those discussed here. Nonetheless, these ex- 
periments demonstrate that microcavities could support 
exotic ordered states that have proved elusive in equilib- 
rium. 



VI. SUMMARY 

We have calculated the linear response of a microcavity 
with a non-equilibrium population of excitons. The pop- 
ulation produces new collective modes, which are analogs 
of the Cooper pairing mode in superconductors. We have 
shown that these modes are visible as peaks in the opti- 
cal spectra. By considering the growth exponents of these 
collective modes we have found a phase diagram for the 
dynamical condensation. In a microcavity with a contin- 
uum of in-plane wavevectors there can be multiple un- 
stable modes of different wavevectors. For some parame- 
ters the dominant (and, for sufficiently negative detuning, 
only) instabilities can occur at a non-zero wavevector. In 
these regimes the microcavity will develop a condensate 
with spatial structure, signaled by coherent emission at 
an angle to the cavity normal. 
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Appendix A: Analytical pump solution 

Ref. 15 gives the results of numerical simulations of 
Eqs. ([l}|3]), driven by a linearly-chirped Gaussian pump 
pulse. These simulations show that there are parame- 
ter regimes in which the dynamics separates into a fast 
pumping stage, followed by a slower condensation stage. 
In this appendix we present an approximate analytical 
solution to Eqs. ([l}|3]) which gives the population and 
(negligible) polarization at the end of the pump pulse. 
This solution forms the starting point for the dynamics 
discussed in the body of the paper. 

The full numerical solutions in Ref. [K] show that the 
only significant polarization during the pumping is at the 
pump wavevector p. Moreover, this polarization can be 
seen to be small compared with the applied pump field 
F. Thus during the pumping we may neglect the second 
term in Eq. (IT| for all wavevectors. With this approxi- 
mation, Eqs. flrpl reduce to an ensemble of independent 
two-level systems, driven by a field ip^ which is the ex- 
ternally applied field F filtered by the cavity response. 
For pumping at high angles, outside the stop band of 
the mirrors, ip'^ is proportional to the pump pulse. Thus 



Eqs. (T][3l become the Bloch equations 




E- A{t) 
-E + A{t) gn{t) 
-g^{t) 



Py 
.Do J 



where we have defined 



-i/* A(t) dt' 



P' - 
P' 



(Al) 

(A2) 
(A3) 
(A4) 



Note that the collective polarizations are at the pump 
wavevector, while the collective inversion is spatially uni- 
form. 

For a model pump pulse of the form 



gQ{t) — — sech -. 

T T 



(A5) 



— tanh 

r T 



Eq. (Al) has an analytical solutiorP^. The form of the 



population Dq at times t ^ t after the pulse is: 



PpjE) 
v{E) 



= 2- 



cosh 



cosh 
\^{{E-uo)r 



1 



COS' 



a)]cosh[f((S-iyo)T + a)] 



- 1. 



(A6) 
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In the limit Qq > a ^ 1 the distribution becomes 



2nF{E - - vo) (1 - npiE - A^- - ^o)) - 1, 

(A7) 

where np{E) is a Fermi distribution with temperature 
l/{kBTTT) and the chemical potentials /i± — i/q zL ^. If 
the density of states i^(-B) is sufficiently small at ener- 
gies below /i_ then this lower edge is irrelevant. The 
occupation function Dq (E) is then equivalent to an equi- 
librium Fermi distribution with /i = fi^. In this pa- 
per we consider parameters where this applies, choosing 
r/p = -30/t, flo = 18 and a = 17.5. 



Since the dynamics during the pumping, Eq. ( Al ), in- 



Pk^p remains zero. For the polarization at the pump 
wavevector, the analytical solution gives a window of en- 
ergies ^ r in which there is a non-zero polarization after 
pumping. However, in the absence of an external field 
and with ipp « 0, as is the case after pumping, the sub- 
sequent evolution of the polarization is free. As a result, 
the total polarization Pp = J Pp{E)dE decays by free 
induction decay, and so may be neglected after a time of 
order r. The numerical work of Ref. IT51 showed that for 
suitable parameters the preparation of this state, includ- 
ing the free induction decay of the remnant polarization, 
finishes before the dynamics discussed in the main body 
of this paper takes place. 



volves only Pp, the polarization at any other wavevector 
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